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Abstract 

A search for squarks in /2-parity violating supersymmetry is performed in e ± p col- 
lisions at HERA using the HI detector. The full data sample taken at a centre-of-mass 
energy yfs = 319 GeV is used for the analysis, corresponding to an integrated luminosity 
of 255 pb" 1 of e + p and 183 pb" 1 of e~p collision data. The resonant production of squarks 
via a Yukawa coupling A' is considered, taking into account direct and indirect ii-parity 
violating decay modes. Final states with jets and leptons are investigated. No evidence for 
squark production is found and mass dependent limits on A' are obtained in the framework 
of the Minimal Supersymmetric Standard Model and in the Minimal Supergravity Model. 
In the considered part of the parameter space, for a Yukawa coupling of electromagnetic 
strength A' = 0.3, squarks of all flavours are excluded up to masses of 275 GeV at 95% 
confidence level, with down-type squarks further excluded up to masses of 290 GeV. 
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1 Introduction 



The ep collider HERA is ideally suited to search for new particles coupling to electrorQ-quark 
pairs. In supersymmetric (SUSY) models with i?-parity violation (Jjtp), squarks can couple 
to electrons and quarks via Yukawa couplings A'. At HERA, squarks with masses up to the 
electron-proton centre-of-mass energy, yfs = 319 GeV, could be produced resonantly via the 
fusion of the incoming 27.6 GeV electron and a quark from the incoming 920 GeV proton. 
Squark decays typically result in a number of high energetic particles in the final state, thus 
several complementary multi-lepton and multi-jet topologies are investigated. The data used in 
this analysis correspond to an integrated luminosity of 255 pb" 1 for e + p collisions and 183 pb -1 
for e~p collisions which represents the full data sample collected at y/s = 319 GeV. For the 
e~p sample, this represents an increase of a factor of thirteen compared to the previous HI 
analysis JT), while for the e + p sample this corresponds to a factor of four. The search presented 
here supersedes the results previously obtained by HI [Q]|2]|. Complementary direct searches for 
Rp SUSY have been carried out at the LEP e + e~ collider [[3HH and at the Tevatron pp collider [[51 
HU. Indirect constraints from low energy precision observables are also available 11714101 . 



2 Phenomenology and Monte Carlo Simulation 

2.1 Production of squarks in H-parity violating supersymmetry 

Supersymmetric extensions of the Standard Model (SM) introduce new elementary particles 
which are the superpartners (sparticles) of SM particles but differ in spin by half a unit. A 
new quantum number R p = (—iy B + L + 2S j s defined, denoted i?-parity, where B is the baryon 
number, L the lepton number and S the spin of a particle. For particles R p = 1 and for their 
supersymmetric partners R p = —1. Most of the collider searches focus on SUSY models 
that conserve i?-parity, allowing only pair-production of sparticles. However, the most gen- 
eral supersymmetric theory that is renormali sable and gauge invariant with respect to the Stan- 
dard Model gauge group does not impose i?-parity conservation. Couplings between two SM 
fermions and a squark (q) or a slepton (/) are then possible, allowing the single production of 
sparticles. The R\, Yukawa couplings responsible for squark production at HERA originate from 
a lepton number violating term X'^ k LiQjDk in the superpotential, where i,j and k are family 
indices. L { , Qj and D k are superfields, which contain the left-handed leptons, the left-handed 
up-type quarks and the right-handed down-type quarks, respectively, together with their SUSY 
partners. Non-vanishing couplings \\- k allow the resonant production of squarks at HERA via 
eq fusion [fTTTl . Feynman diagrams of these processes are shown in figured] The values of 
the couplings are not fixed by the theory but are required to be small to conform with present 
observations. For simplicity, it is assumed here that one of the A' ljjt couplings dominates over 
all the other trilinear ft p couplings. At high Bjorken-x the density of antiquarks in the proton is 
significantly smaller than that of the valence quarks. Hence e~p scattering gives sensitivity to 
the couplings A' llfe (k = 1, 2, 3) which dominate the production of d^-type squarks (i.e. the su- 
perpartners dft, sr and of down-type quarks). The dominant contribution to the production 

'In the following the generic term electron refers to both electron and positron unless explicitly stated otherwise. 
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cross section is thus approximately proportional to Xf lk ■ u(x) where u(x) gives the probability 
to find a u quark in the proton carrying the momentum fraction x = M~ / s, where M| is the 
squared mass of the produced squark. By contrast, e + p scattering provides sensitivity to the 
couplings (j = 1, 2, 3) which dominate the production of w^-type squarks (i.e. the super- 
partners ul, cl and ti of up-type quarks). Here the dominant contribution to the production 
cross section is approximately proportional to Xf^ ■ d(x). Due to the larger u quark density in 
the proton at large x with respect to the d quark density, larger production cross sections are 
expected in e~p interactions for identical couplings and squark masses. 

Signal cross sections are obtained in the narrow width approximation by using the leading 
order amplitudes given in lfi~2ll . corrected to account for next-to-leading order QCD effects 
using multiplicative correction factors lfP3ll . The parton densities are evaluated at the hard scale 
Mf. 

2.2 Final states from squark decays 

In |? p SUSY all sparticles are unstable. Squarks can decay directly via the Yukawa coupling A' 
into SM fermions. The c^-type (k = 1,2,3) squarks can decay via the coupling A' llfc either into 
e~ + uox v e + d, while the u^-type (j = 1, 2, 3) squarks decay via the coupling X[ - x into e + + d 
only, as illustrated in figure [TJ Squarks may also decay via R p conserving gauge couplings as 
illustrated in figure [2] The u^-type squarks can undergo a gauge decay into states involving 
a neutralino x°i (i — 1> 2, 3, 4), a chargino xt = h 2) or a gluino g. In contrast, d^-type 
squarks mainly decay to x°i or 9 an d decays into charginos are suppressed IfTTTl . 

The squark decay chains analysed in this paper are classified by event topology [UHl. This 
classification relies on the number of isolated electrons, muons and hadronic jets in the final 
state, and on the presence of missing energy (indicating undetected neutrinos). The channels 
labelled eq and vq are the squark decay modes that proceed directly via|? p couplings resulting in 
event topologies with an isolated electron or neutrino and a single jet. The remaining channels 
result from the gauge decays of the squark and are characterised by multijet (M7) final states 
with additional leptons. The channels labelled eMJ and uMJ involve one or two gauginos (x or 
g) in the decay cascade. In the eMJ channel e + and e~ are possible in the final state, such that 
with respect to the incident lepton charge, a "right" (same sign) charge eMJ(RC) and a "wrong" 
(opposite sign) charge eMJ(WC) channel are distinguished. Channels with an electron and an 
additional charged lepton £ (where i = e,/x) denoted by eeMJ, e/jMJ or neutrinos and a 
charged lepton evMJ, ufxMJ (generically written as e£MJ and vlMJ) necessarily involve two 
gauginos. Decay patterns involving more than two gauginos are kinematically suppressed and 
are therefore not explicitly studied here. Processes leading to final states with tau leptons are 
not expected to increase the sensitivity of this analysis and are not explicitly investigated. A 
dedicated search for isolated tau leptons in HI shows good agreement with the SM lfi~4~l . 

2.3 Event simulation 

For each of the signal topologies described above a dedicated Monte Carlo (MC) simulation is 
done. For the direct lepton-quark decay channels eq and vq, as shown in figure Q3 events are 
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generated using LEGO lTT5i For the gauge decays of squarks (figure© events are generated 
using SUSYGEN3 |[T6t 

To allow a model independent interpretation of the results, the squark decay processes are 
simulated for a wide range of masses of the sparticles involved. The final states contain only 
SM fermions (/) considered as massless. The squark mass is varied from 100 GeV to 290 GeV 
in steps of typically 25 GeV. For gauge decays of squarks involving a gaugino, which decays 
directly via |? p , the process q — > qx\ is generated for Xi masses ranging between 30 GeV 
and Mq. In order to study cascade gauge decays which involve two gauginos, the processes 
<? — ► QXi ~~ ^ QXiff' an d q qx® QXiff' are generated for xf an d X2 masses ranging 
between 40 GeV and M„, and for v? masses between 30 GeV and NL + or M v a. The masses of 
the gauginos are varied in steps of approximately 10 GeV. The lower mass values for squarks 
and gauginos are motivated by the exclusion domains resulting from fjL p SUSY searches at 
LEP l3l|4l. The mass intervals are sufficiently small to allow linear interpolation of signal 
detection efficiencies as a function of the masses of the sparticles involved. 

SM processes may mimic the characteristics of the final states of squark decays. This SM 
background is dominated by neutral current (NC) and charged current (CC) deep inelastic scat- 
tering (DIS), with additional small contributions from photoproduction, single W boson pro- 
duction and lepton pair production. The RAPGAP IfTTl event generator, which implements the 
Born level, QCD Compton and boson-gluon fusion matrix elements, is used to model inclusive 
NC DIS events. The QED radiative effects arising from real photon emission from both the 
incoming and outgoing electrons are simulated using the HERACLES [18J program. Direct 
and resolved photoproduction of jets and prompt photon production are simulated using the 
PYTHIA |fT9ll event generator. The simulation is based on Born level scattering matrix ele- 
ments with radiative QED corrections. In RAPGAP and PYTHIA, jet production from higher 
order QCD radiation is simulated using leading logarithmic parton showers and hadronisation is 
modelled with Lund string fragmentation ll2~0~l . Inclusive CC DIS events are simulated using the 
DJANGO ETi program, which includes first order leptonic QED radiative corrections based on 
HERACLES. The production of two or more jets in DJANGO is accounted for using the colour 
dipole model (221 • The leading order MC prediction of processes with two or more high trans- 
verse momentum jets in NC DIS, CC DIS and photoproduction is scaled by a factor of 1.2 to 
account for the incomplete description of higher orders in the MC generators [|23ll24l . Contri- 
butions arising from the production of single W bosons and multi-lepton events are modelled 
using the EPVEC E51 and GRAPE ll26ll event generators, respectively. 

Generated events are passed through a GEANT 11271 based simulation of the HI apparatus, 
which takes into account the actual running conditions of the data taking. Simulated events are 
reconstructed and analysed using the same program chain as is used for the data. 

3 Experimental Method 
3.1 HI detector 

A detailed description of the HI experiment can be found elsewhere ll28l . Only the detec- 
tor components relevant to this analysis are briefly described here. A right-handed cartesian 
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coordinate system is used with the origin at the nominal primary ep interaction vertex. The 
proton beam direction defines the positive z axis (forward direction). The polar angle 6 and 
the transverse momenta P T of all particles are defined with respect to this axis. The azimuthal 
angle <p defines the particle direction in the transverse plane. The pseudorapidity is defined as 
f] = — In tan |. 

The Liquid Argon (LAr) calorimeter ||2~9| covers the polar angle range 4° < 6 < 154° with 
full azimuthal acceptance. The energies of electromagnetic showers are measured in the LAr 
with a precision of a(E)/E ~ 11%/y/E/GeV® 1% and hadronic energy depositions with 
a(E) /E ~ 50%/ ^E /GeV © 2%, as determined in test beam measurements J30l[3ll- A lead- 
scintillating fibre calorimeter (SpaCal) |[32ll covering the backward region 153° < 6 < 178° 
completes the measurement of charged and neutral particles. For electrons a relative energy 
resolution of a(E)/E ~ 7%/ \/E/ GeV © 1% is reached, as determined in test beam measure- 
ments J33]|. The central (20° < 9 < 160°) and forward (7° < 9 < 25°) inner tracking detectors 
are used to measure charged particle trajectories and to reconstruct the interaction vertex. The 
LAr calorimeter and inner tracking detectors are enclosed in a super-conducting magnetic coil 
with a field strength of 1.16 T. From the curvature of charged particle trajectories in the mag- 
netic field, the central tracking system provides transverse momentum measurements with a 
resolution of o Pt /P t = 0.005P r /GeV © 0.015 11341 . The return yoke of the magnetic coil is 
the outermost part of the detector and is equipped with streamer tubes forming the central muon 
detector (4° < 9 < 171°). In the very forward region of the detector (3° < 9 < 17°) a set of 
drift chambers detects muons and measures their momenta using an iron toroidal magnet. The 
luminosity is determined from the rate of the Bethe-Heitler process ep — > e/ry, measured using 
a photon detector located close to the beam pipe at z = — 103 m, in the backward direction. 

3.2 Particle identification and event reconstruction 

Electromagnetic particle (electron and photon) candidates are identified as compact and isolated 
clusters of energy in the electromagnetic part of the LAr calorimeter. Electron candidates are 
identified as electromagnetic particle candidates with an associated track. Identification of muon 
candidates is based on a track in the inner tracking detectors, associated to a signal in the muon 
system. Tracks and calorimeter deposits not identified as originating from isolated electro- 
magnetic particles or muons are combined into cluster-track objects to reconstruct the hadronic 
final state 11351 . Jets are reconstructed from these objects using an inclusive k T algorithm B6II37I 
with a minimum P T of 4 GeV. The missing transverse momentum P™ 1SS , which may indicate the 
presence of neutrinos in the final state, is derived from all reconstructed particles in the event. A 
neutrino four-vector P u is reconstructed by exploiting momentum and energy conservation. The 
transverse momentum of the neutrino is reconstructed by assuming one neutrino with significant 
energy in the event P% = P™ 1SS . The energy of the neutrino is then reconstructed exploiting the 
energy and longitudinal momentum balance: (E l — P*) + [E v — P%) = 2P° = 55.2 GeV, 
where the sum runs over all detected particles, P z is the momentum along the proton beam axis 
and P° denotes the energy of the incident electron. 

For further selection the following observables are used which in SM DIS events correspond 
to the Lorentz-invariant quantities: inelasticity y, negative four-momentum transfer squared Q 2 
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and Bjorken's scaling variable x. They can be reconstructed as: 

_ _ P e (l-cosfl e ) 2 _ _?fe_ _ 01 

where the polar angle 9 e , energy E e and transverse momentum Pf; of the electron with the 
highest Pt found in the event are used. If no electron is reconstructed in the event similar 
quantities can be calculated using the Jacquet-Blondel method [38] from the hadronic final 
state variables: 

yh — o Z7T) ' ~ i — i 

where Pr,h is the transverse momentum of the hadronic final system calculated from the before 
mentioned cluster-track objects. The sum {E — P z )h mns over energies and momenta of 
jets only. 



3.3 Trigger and data quality 

All data events used for this search are triggered by the LAr calorimeter ||39l . Events with 
an electromagnetic deposit in the LAr with an energy greater than 10 GeV are detected with 
an efficiency close to 100% BD1 . Events are also triggered by hadronic jets, with a trigger 
efficiency above 95% for a jet transverse momentum P^ et > 20 GeV and almost 100% for 
> 25 GeV iHTTl . For events with missing transverse energy of 20 GeV, the trigger efficiency 
is about 90% and increases to above 95% for missing transverse energy above 30 GeV ll4~2ll . 

In order to remove background events induced by cosmic showers and other non-ep sources, 
the event vertex is required to be within 35 cm in z of the mean position for ep collisions. In 
addition, topological filters and timing vetoes are applied ll43l in order to reject beam related 
and cosmic background. 



4 Data Analysis 

The event selection is carried out in several exclusive analysis channels. The resulting event 
rates for each channel and the range of efficiencies for the selection of signal events are given 
in tabled! 



4.1 Electron-jet final state eq 

The final state of a squark decaying into an electron and a high P T jet is identical to the 
NC DIS signature at high x and Q 2 . For the signal the reconstructed invariant mass dis- 
tribution M e = ^/x e s shows a resonance peak at the nominal squark mass with a resolution 
SM e = 4 — 10 GeV depending on the mass of the squark. Differences in the M e and y e distri- 
butions of the two processes allow to discriminate them statistically. Squarks produced in the 
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s-channel decay isotropically leading to a flat da/dy distribution, whereas for NC DIS events 
a distribution proportional to 1/y 2 is expected. 

Events are selected by requiring P^ iss < 15 GeV and 40 GeV < - P z ) < 70 GeV 
where the sum runs over all particles in the final state. An isolated electron with Pf. > 16 GeV 
in the region 5° < 9 e < 145°, Q\ > 2500 GeV 2 and y e < 0.9 is required. The high y region is 
excluded to reduce background contributions arising from photoproduction events. An M e de- 
pendent cut on y e is determined by minimising the expected limit from signal and SM MC. The 
cut ranges from y e > 0.5 for masses around 100 GeV to y e > 0.2 for masses around 290 GeV. 
In order to remain exclusive with respect to other channels, events with muon candidates with 
P% > 5 GeV or two jets with Pf * > 15 GeV are rejected. 

The M e spectra after this selection for data and SM background are shown in figures [3k 
and|4k for e~p data and e + p data, respectively. For both data samples no significant deviation 
from the SM expectation is observed. In the e~p sample, 3121 events are observed while the 
SM expectation yields 3215 ± 336. In the e + p data sample 2946 candidate events are found 
compared to 2899 ± 302 expected from SM processes. 

4.2 Neutrino-jet final state vq 

Squarks decaying into a neutrino and a high P T jet lead to the same signature as CC DIS events 
with high missing transverse momentum. Similarly to the eq channel, the resonant s-channel 
production and isotropic decay allows a statistical separation of signal and background. The 
resolution SMh of the reconstructed invariant mass Mh = yjxh~s varies between 12 and 22 GeV 
depending on the mass of the squark considered. 

The presence of a neutrino in the event is required by imposing P™ 1SS > 30 GeV and 
52(E- P z ) < 50 GeV. The phase space is restricted to Q\ > 2500 GeV 2 and y h < 0.9. 
Similar to the eq channel, a cut on y h dependent on the reconstructed mass M h is applied. The 
cut ranges from y^ > 0.3 for masses around 100 GeV to yh > 0.1 for masses around 290 GeV. 
In order to remain exclusive with respect to other channels, events with any electron or muon 
candidate with P^' M > 5 GeV or events containing two jets with Pf > 15 GeV are rejected. 

Only d,R- type squarks, which are produced dominantly in e~p collisions, can undergo direct 
decay leading to a vq final state. For selected events the M h spectrum of this data set and of the 
simulation of SM background events is shown in figure (3h- No significant deviation from the 
SM expectation is found. In data 2858 events are observed while 2983 ± 358 are expected from 
SM processes. 

4.3 Electron-multijet and electron-lepton-multijet final states 

Common preselection for eMJ(RC) , eMJ(WC) , eeMJ, e^iMJ and euMJ 

Squarks decaying via neutralinos or charginos are expected to have a higher multiplicity of jets 
and leptons in the final state. The signatures correspond to final states detectable in higher order 
NC DIS processes. However, as heavy particles are boosted forward, in events with squarks the 



9 



decay products are mainly emitted into the forward part of the detector. This feature is used to 
distinguish between the signal and SM background. 

A common preselection is applied for the eMJ(RC), eMJ(WC), eeMJ, efiMJ and evMJ 
channels: At least one isolated electron with Fj, > 6 GeV and E e > 11 GeV in the region 
5° < 9 e < 110° is required. In addition, the condition y e > 0.3 is used to reduce the back- 
ground from NC DIS. Central electrons with 30° < 9 e < 110° are required to have a well mea- 
sured track associated to the cluster and the distance of closest approach between the track 
impact point and the centre-of-gravity of the cluster should not exceed 12 cm. Furthermore 
the energy of the track is required to match the energy of the associated cluster according 
to Cluster /-^track > 0-5- At least two jets with P^ ct > 15 GeV in the jet polar angle range 
7° < 6jct < 145° are also required. 

By requiring Q 2 > 1000 GeV 2 the steep decrease of the NC DIS cross section with increas- 
ing Q 2 is exploited. This corresponds to an implicit upper cut on the electron polar angle as 
Q 2 is strongly correlated with the polar angle of the scattered electron. For signal events at 
least one high Pt particle is expected to be emitted in the forward direction, therefore the high- 
est Pt electron or one of the two highest P T jets has to fulfill 6> e j et < 40°. Moreover, of the 
two highest Pt jets, the one with the largest polar angle, #backw> must satisfy the condition 
#backw < (He — 0.3) • 180°, separating efficiently signal events from NC DIS background fl4l . 

For selected events an invariant mass M roc is calculated as M rec = y/ 4B°Q^ E{ — E®), 
where the energies Ei of electrons, muons and jets found in the event with P^ et > 5 GeV are 
included in the sum. The resolution SM rec of this method ranges between 6 and 10 GeV de- 
pending on the mass of the squark considered. 

Electron-multijet final state eMJ(RC), eMJ(WC) 

Decays of squarks via gauginos are likely to produce a single isolated electron and multiple 
jets in the final state. Events produced in e ± p collisions with exactly one electron and multiple 
jets in the final state are also expected in the SM, where in general the measured charge of the 
electron corresponds to that of the incident electron. A selection channel labelled "right" (same 
sign) charge eMJ (RC) is used for events fulfilling this criterion. A selection channel where 
the electron charge is identified as opposite to the incident electron, denoted "wrong" (opposite 
sign) charge eMJ (WC), represents therefore a powerful test of the SM and is expected to be 
essentially background free. The distinction between RC and WC eMJ events is based on the 
curvature of the electron track measured in the central tracking system. Events are allocated to 
the WC channel if the electron is found in the central region 30° < 9 e < 110° and its charge is 
measured to be opposite to that of the incident electron, with a charge significance greater than 
two standard deviations [44]. Otherwise the event is assigned to the RC channel. 

Events are selected from the common preselection described above by requiring in addition 
pmiss < Q e y an( j 4Q Q e y < Y^(E - P z ) < 70 GeV since no neutrinos are expected in these 
channels. To ensure that the selection channels are exclusive, no additional electron or muon 
candidate with P T ' ^ > 5 GeV may be present in the event. 

In the eMJ (WC) channel no event is observed in the e~p data while 1.3 ± 0.3 are expected 
from SM processes and one candidate event is observed in the e + p data while 0.6 ± 0.4 are 
expected. For e + p collisions the mass spectrum is shown in figure Hb- 
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For events in the eMJ (RC) channel an M rec dependent cut on y e is applied to increase 
the sensitivity for signal events. The cut ranges from y e > 0.7 for masses around 100 GeV to 
y e > 0.5 for masses around 290 GeV. The M rec distributions for data and simulation are shown 
in figures |3fc and0t for e~p and e + p collisions, respectively. No significant deviation from the 
SM expectation is observed. In total 147 events are observed in the e~p data while the SM 
simulation yields 158.3 ± 23.9 and in the e + p data 140 events are observed for an expectation 
of 146.0 ± 21.4 from SM processes. 

Electron-lepton-multijet final states, eeMJ, e^iMJ, euMJ 

Final states from squark decays may contain more than one isolated lepton if the decays pro- 
ceed via cascades of gauginos. In addition to the common preselection the eeMJ and e/xMJ 
channels require either an additional electron with the same criteria as described in the common 
preselection or an isolated muon with P£ > 5 GeV in the polar angle range 10° < 9^ < 110°. 
After applying this selection the SM background expectation is very low in these channels. For 
the e/xMJ channel there are no candidate events observed in the e~p and e + p collision data for 
a SM expectation of 0.03 ± 0.02 and 0.03 ± 0.03, respectively. In the e~p collision data no can- 
didate event for a SM expectation of 1.5 ± 0.5 is observed in the eeMJ channel and two events 
in the e + p data are observed compared to a SM expectation of 1.7 ± 0.5. The mass spectrum 
for the eeMJ channel in the e + p data is shown in figure HJl. 

In the evMJ channel a neutrino is expected in the final state, therefore in addition to the 
common preselection, large missing transverse momentum P™ 1SS > 15 GeV is required. Due to 
the presence of the neutrino ^2(E — P z ) is significantly reduced causing to be substantially 
smaller than y e , while in NC DIS events y^ « y e is expected. Thus a cut y e (y e — Uh) > 0.04 is 
used to discriminate the SUSY signal from background events ll44l . Exclusivity with respect to 
the eeMJ and efiMJ channels is achieved by rejecting events containing an additional electron 
or muon with > 5 GeV. The method used to reconstruct a mass M reCjl/ for selected events 
taking the energy of the neutrino into account is explained in section |4~4| 

In the e~p collision data three events are observed in the evMJ channel while 5.6 ±1.2 are 
expected and in e + p collision data five events are observed while 8.2 ± 2.0 are expected from 
SM processes. The mass spectra are shown in figures [3]1 and|4^ for the e~p data and e + p data, 
respectively. 

4.4 Neutrino-multijet and neutrino-muon-multijet final states 
Common preselection for uMJ and v^iMJ 

Squark decays with single or multiple neutrinos produced via neutralino or chargino decays can 
result in final states similar to that of higher order CC DIS processes. 

A substantial missing transverse momentum P™ 1SS > 26 GeV is required and at least two jets 
must be found with P^ > 15 GeV in the range 7° < 6 ]Ct < 145°. No electron candidate with 

> 5 GeV is allowed. A cut J2(E — P z ) < 50 GeV is used to ensure the neutrino energy is 
positive. 
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For each selected event a squark mass M reC)J/ is calculated as M rec u = \/ii^Q^ Ei — Efj 
where the sum includes the energies of the reconstructed neutrino, electrons, muons and jets 
with P^ ct > 5 GeV in the event. This method assumes that all missing energy is carried by a 
single neutrino and yields a resolution 5M reCjU of about 15 to 20 GeV depending on the squark 
mass. 

Neutrino-multijet final state uMJ 

Squark decays via gauginos are likely to produce final states with multiple jets and a single 
neutrino. Events are selected in the uMJ channel if no muon candidate is found. 

A cut on y h dependent on the reconstructed mass M reCjI/ is applied to enhance the signal. The 
cut ranges from y h > 0.5 for masses around 100 GeV to y h > 0.4 for masses around 290 GeV. 
The M rec „ spectra are shown in figures [3k and|4j for e~p and e + p collision data and SM back- 
ground simulation. In the e~p data 204 candidate events are selected while 235.5 ± 63.3 are 
expected and in the e + p data 113 candidate events are selected while 134.0 ± 33.8 are expected 
from SM processes. 

Neutrino-muon-multijet final state ufiMJ 

If an isolated muon with > 5 GeV in the polar angle range 10° < 9^ < 110° is found in an 
event in the common vMJ selection, the event is classified as ujiMJ candidate. No candidate 
events are found, in the e~p or in the e + p collision data, in agreement with the SM expectations 
of 0.04 ± 0.02 and 0.06 ± 0.03, respectively. 

4.5 Systematic uncertainties 

The following experimental systematic uncertainties are considered: 

• The uncertainty on the electromagnetic energy scale varies depending on the polar an- 
gle from 0.7% in the central region to 2% in the forward region []42|. The polar angle 
measurement uncertainty of electromagnetic clusters is 3 mrad. 

• The jet energy scale is known within 2% ll42l . The uncertainty on the jet polar angle 
determination is 10 mrad. 

• The luminosity measurement has an uncertainty of 3%. 

The effects of the experimental systematic uncertainties on the SM expectation and signal effi- 
ciencies are determined by varying the corresponding experimental quantities within one stan- 
dard deviation in the MC samples and propagating the variations to the final distributions. The 
resulting experimental uncertainties are determined for each analysis channel individually and 
added in quadrature. In the eq channel the uncertainty on the overall SM event yield was found 
to be 3%, while in the vq channel an uncertainty of 7% is determined. In the eMJ (RC and WC) 
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channels the resulting systematic uncertainty amounts to 4%, in the efiMJ and evMJ channels 
to 7%, in the eeMJ channel to 10% and in the uMJ and vfiMJ channels to 20%. 

Additional model uncertainties are attributed to the SM MC event generators described in 
section[23l A conservative error of 10% is attributed to NC (RAPGAP) and CC (DJANGO) DIS 
processes with only one high P T jet. To account for the uncertainty on higher order QCD correc- 
tions, an uncertainty of 15% is attributed to NC DIS and photoproduction processes (PYTHIA) 
with at least two high P T jets. The normalisation uncertainty of CC DIS processes with at least 
two high P T jets is estimated to be 20% lf42~l . A 5% uncertainty is attributed to the contribu- 
tion from multi-lepton events (GRAPE) and a 15% uncertainty on the production of single W 
bosons (EPVEC). These uncertainties include contributions from the proton parton distribution 
functions and from missing higher order QCD corrections. The total error on the SM prediction 
is determined by adding the effects of all model and experimental systematic uncertainties in 
quadrature. 

For the signal cross section further uncertainties arise from the determination of signal effi- 
ciencies (10% due to available MC statistics), the theoretical uncertainty on the squark produc- 
tion cross section (7% for low squark masses, up to 50% for the highest masses from the PDF 
uncertainty) and an uncertainty due to the scale at which the PDFs are evaluated (7%) Il44l . 



5 Exclusion Limits 

No significant deviation from the SM expectation is observed in any channel. Consequently the 
observations in all analysis channels are combined to set constraints on various supersymmet- 
ric models. Exclusion limits are obtained on the production of squarks parameterised by the 
strength of the ]jt p couplings A' ljX and A' llfe and dependent on the mass of the squark. 

5.1 Procedure 

For the interpretation of the results a version of the Minimal Supers ymmetric Standard Model 
(MSSM) is considered where the masses of the neutralinos, charginos and gluinos are deter- 
mined via the usual parameters: the "Higgs-mass" term fi, which mixes the Higgs superfields; 
the SUSY soft-breaking mass parameter M 2 ; and the ratio of the vacuum expectation values of 
the two neutral scalar Higgs fields tan/3 [|45l . The parameters are defined at the electroweak 
scale. 

A set of parameters (tan f3, fi, M 2 ) together with the sfermion masses and a coupling \' ljk 
define a supersymmetric scenario where the masses of the gauginos and the branching ratios 
for squark decays into the different final state topologies are fixed and can be obtained using 
the SUSYGEN3 [16] package. The branching ratios for the specific parameters of the model 
are taken into account in the combination ll46l . A sliding mass window technique is used in 
channels with high contributions of irreducible SM background (eg, uq, eMJ and uMJ) to 
improve the signal to background ratio for the squark mass examined. The width of the mass 
window is determined by minimising the expected limit, and increases towards high squark 



13 



masses, reflecting the corresponding mass reconstruction resolution. The small efficiency losses 
due to the finite mass window width are taken into account. A 95% confidence level (CL) upper 
limit cr lim on the squark production cross section compatible with the simultaneous observation 
in all channels is derived using a modified frequentist approach based on Likelihood ratios ll46ll . 
Sets of model parameters leading to signal cross sections above ay m are excluded. 

If the squark width is non-negligible, in particular for squark masses approaching the kine- 
matic limit, the production cross section decreases at the resonance peak and contributions from 
the lower tail of the squark mass distribution become important, enhanced by the rapid increase 
of proton parton distributions at low Bjorken-x ll44l . This is taken into account by generating 
events for negligible squark widths to determine signal efficiencies at all masses. The selection 
efficiencies are then corrected for the actual squark width by reducing the efficiency for the 
signal selection accordingly 11441 . 

In the special case of stop and sbottom squark production, namely via|? p couplings A' 131 and 
A' 113 , the mixing of the weak eigenstates t L , t R (and b L , b R ) to the mass eigenstates t\, t 2 (&i, b 2 ) 
via an angle 9$ (#5) becomes important for the calculation of branching ratios and production 
cross sections. The gauge decay via a top quark would lead to decay products different from 
the first two generation squarks, for which the efficiencies are determined. Final states with top 
quarks are not considered explicitly. Since top signals would in any case be present in one of 
the selection topologies, this approach is conservative. 

5.2 Constraints on a phenomenological MSSM 

Constraints are set in a scenario of a phenomenological MSSM P51 where the lightest super- 
symmetric particle (LSP) is the neutralino Xi- Slepton masses Mf are fixed at 90 GeV, close 
to the lowest mass bound from fji v sfermion searches at LEP [3j and squark masses are treated 
as free parameters. For higher slepton masses only very small degradations in the derived con- 
straints are expected JT). 

For a single point in the parameter space, characterised by // = —200 GeV, M 2 = 80 GeV 
and tan (3 = 2, constraints on the strength of the couplings depending on the mass of the 
squark are derived for d R (k = 1,2) (figure [5^) and u 3 L (j = 1,2) (figure O?) production. 
The HERA sensitivity allows tests of A' values as low as 10~ 2 for squark masses of 100 GeV. 
For high squark masses the sensitivity degrades since the production cross section decreases 
strongly. The limits from the previous HI analysis on a smaller data sample HI are also indi- 
cated. 

This choice of parameters leads to a dominant photino (7) component to the neutralino's 
composition. As a consequence, gauge decays are likely to result in charged leptons in the 
final state. The branching ratios into the decay topologies are shown at the observed limit 
for d\(k = 1,2) (figure Et) and u j L (j = 1,2) (figure [5)3) production. For d\{k = 1,2) 
production the channels eMJ (RC) and eMJ (WC) each contribute about 40% over a wide range 
of squark masses and only 10% of squark decays appear in the vMJ channel. For squark masses 
approaching the kinematic limit of the centre-of-mass energy the lepton-quark channels eq and 
vq begin to dominate the decays of the squarks, because gauge decay modes become negligible 
at high values of the couplings. Over the whole mass range the sum of analysed branching 
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ratios is close to 100%. For u 3 L (j = 1,2) production the channels e£MJ have the highest 
branching ratio over a wide mass range at the observed limit. 

A different point in the parameter space yields a complementary scenario with the choice 
of // = 200 GeV, M 2 = 150 GeV and tan (3 = 2. Again the neutralino x°i is the LSP but its 
composition is now dominated by a zino (Z) component. Squark decays are now more likely to 
produce neutrinos in the final state. The constraints on the couplings depending on the mass of 
the squark are shown for d R (k = 1,2) (figured) and u J L (j = 1, 2) (figure^) production and 
are of the same order of magnitude as in the photino scenario. Branching ratios at the observed 
limits show dominant contributions from the uMJ and vlMJ channels (figures |5fc and[6]l). 

These two scenarios illustrate the sensitivity for various model configurations achieved by 
the combination of the complementary search topologies. The sensitivity of the analysis is 
explored in a scan of the MSSM parameters. The parameters M 2 and fi are varied in the range 
70 GeV < M 2 < 350 GeV and -300 GeV < fi < 300 GeV for tan/3 = 6. Parameter 
sets leading to a scalar LSP or to LSP masses below 30 GeV are not considered. The latter 
restriction, as well as the lower boundary of the M 2 range, are motivated by the exclusion 
domains resulting from gaugino searches in $ p SUSY at LEP flU. Figures [7] and [8] show the 
resulting constraints on the couplings as a function of the squark mass. The region of values 
excluded at 95% CL for the couplings in all scenarios and the best exclusion limit achieved 
in all scenarios are indicated for first and second generation squarks d R , sr (figure [7^) and 
ul, &l (figure [8k) as well as for third generation squarks bjt (figure [7b) and II (figure [8b). 
The resulting exclusion domains are compared to the previous HI results H). Constraints on 
the ftp couplings are also available as indirect limits from low energy experiments probing 
virtual squark contributions [10]. The production of up-type and down-type squarks via the 
A' m coupling is strongly constrained by the non-observation of neutrinoless double beta decay 
(J3f30v) [0[TO]|. The best indirect limit on the couplings A' 112 and A' 113 results from tests of 
charged current universality (CCU) flUEO! and can be compared to the direct limits obtained in 
this analysis for \' llk in figure [7J The best indirect limit on the couplings A' 121 and A' 131 comes 
from atomic parity violation (APV) measurements [HQjl and can be compared to the direct 
limits obtained for A' 1X in figure [8j 

In the part of the parameter space considered here, Yukawa couplings of electromagnetic 
strength A'^ or A' llfc = \ZAna cm = 0.3, are excluded up to masses of 275 GeV at 95% CL for 
up-type squarks and up to masses of 290 GeV for down-type squarks. 

5.3 Constraints on the Minimal Supergravity Model 

Constraints are also obtained on the Minimal Supergravity Model (mSUGRA) PTTl which is a 
complete SUSY model using the assumption of gauge coupling unification and radiative elec- 
troweak symmetry breaking (REWSB) with the choice of 5 parameters: the common mass of 
scalar sparticles m ; the common mass of fermionic sparticles mm; the common trilinear cou- 
pling A ; the ratio of Higgs vacuum expectation values tan (3; and the sign of the Higgs mixing 
parameter ji. The masses of squarks, sleptons and gauginos as well as the branching ratios in the 
analysis channels are determined by the set (m , mi/ 2 , tan/3, sign(jj), A ) for given values of 
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the couplings A' llfe and A' yi . The program SUSPECT 2.1 B3 is used to obtain the REWSB so- 
lution for |/4 1 and to calculate the full supersymmetric mass spectrum. A enters only marginally 
in the interpretation and is set to zero. The parameter ji is taken with negative sign. 

Figures [9] and [TO] show constraints in the m , nii/ 2 plane when values of the couplings are 
assumed to be of the electromagnetic coupling strength A' 1U , = 0.3 or \' x - x = 0.3 for dif- 
ferent values of tan (3. The excluded region typically covers masses of m(u) = 275 GeV, 
mil) = 270 GeV and m(d) = 280 GeV, as indicated in the figures. Complementary constraints 
are obtained by the L3 experiment [4] at LEP and the D0 experiment [6] at the Tevatron which 
exploit di-electron events. The LEP and Tevatron limits are independent of the Yukawa cou- 
pling. For tan (3 = 2, the parameter space is more strongly constrained by the searches for 
gauginos and sleptons at the L3 experiment at LEP, as shown in figures [9] and [jXJ This is the 
only tan (3 value considered in the L3 analysis, although results for higher values are expected 
to be similar Q. Compared to the D0 experiment, the HI limits are more stringent only for low 
values of m for tan (3 = 2, whereas for tan (3 = 6 the domain excluded by HI is considerably 
larger. 

The exclusion limits in figure [9] are very similar for all three flavours of down-type squarks. 
Significant differences are observed between the first two and the third generation of up-type 
squarks. The stronger limit for stop squark production results from strong mixing effects that oc- 
cur for third generation squarks with increasing tan (3, leading to masses for stop squarks lower 
than for first and second generation up-type squarks. The tan (3 dependence of the mSUGRA 
exclusion limits is studied assuming a unified common mass M = m = mi/ 2 and cou- 
plings of electromagnetic coupling strength. This is illustrated in figure [TTJ for the individual 
flavours. For the first two generations of up-type and down-type squarks no dependence on 
tan (3 is observed and values of M < 105 GeV and M < 110 GeV, respectively, are excluded 
over the whole range. In the case of stop squark production, significantly higher values (up to 
M < 148 GeV) are excluded due to the presence of a light stop squark state. This effect is 
also observed for sbottom production, where increased values of tan (3 allow higher values of 
M to be excluded. While there is a steep increase for the limit from stop squark production at 
small tan (3 and a flat plateau over the remaining tan (3 range, the limit from sbottom production 
increases steadily over the range of tan (3 values. The sharp edge in the stop exclusion curve at 
tan (3 ~ 38 for stop production follows from mixing effects in the f sector at high tan (3 leading 
to scenarios with strong contributions from events with r leptons in the final state, which are 
not explicitely considered in the analysis. 



6 Summary 

A search for i?-parity violating production of squarks in 255 pb" 1 of e + p and 183 pb" 1 of e~p 
collisions at HERA is presented. No significant deviation from the Standard Model is observed 
in the study of final state topologies which may result from direct or indirect squark decays. 
Mass dependent limits on the $ p couplings A'^ and \' llk (j, k = 1, 2, 3) are derived within a 
phenomenological version of the MSSM. The existence of w^-type and rf^-type squarks of all 
three generations with masses up to 275 GeV and 290 GeV, respectively, is excluded at the 95% 
CL for Yukawa couplings of electromagnetic strength. These mass limits set the most stringent 
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direct bounds on x and \' llk . For lower squark masses, the results improve the indirect bounds 
set by low-energy experiments. Exclusion limits are also derived in the mSUGRA model, and 
are competitive with and complementary to those derived at the LEP and Tevatron colliders. 
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Table 1 : Total numbers of selected events, SM expectations and ranges of signal efficiencies 
for the squark decay channels considered in e~p and in e + p collisions. The range of signal 
efficiencies gives the extreme values for squark masses ranging from 100 GeV to 290 GeV and 
gaugino masses ranging from 30 GeV up to the squark mass. 

The vq channel is not relevant for e + p data since the w^-type squarks produced in e + p do not 
undergo this decay. Only d R -type squarks, which are produced dominantly in e~p collisions, 
can undergo direct decay leading to a vq final state. The total error on the SM prediction is 
determined by adding the effects of all model and experimental systematic uncertainties in 
quadrature. 
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Figure 1: Feynman diagrams for the single resonant s-channel production of right-handed 
down-type squarks in e~p collisions (a) and left-handed up-type squarks in e + p collisions (b) 
with subsequent decays into SM particles via Yukawa couplings X' llk or A' ljl5 respectively. The 
right-handed down-type squarks can decay either into e~ + u or v e + d, while the left-handed 
up-type squarks decay into e + + d only. 




(a) (b) 

Figure 2: Feynman diagrams for squark decays proceeding via gauginos in the case of right- 
handed down-type squarks (a) and left-handed up-type squarks (b) with subsequent ]jL p decay 
into SM fermions via Yukawa couplings X' llk or A' ljl? respectively. The resulting final states 
may contain multi-leptons and multi-jets. The d R -type squarks decay to Xt d — 1> 2, 3, 4) or 
g, decays into charginos are suppressed while u L -type squarks couple also to charginos xt 
(z = l,2). 
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Figure 3: Reconstructed invariant mass distributions in all selection channels with data (points) 
events from 183 pb -1 of e~p collisions compared to SM MC predictions. The method used for 
the reconstruction (M e , M h , M rec , M rec u ) depends on the analysis channel. The error band gives 
all model and experimental systematic uncertainties on the SM prediction (solid histogram) 
added in quadrature. Error bars of data points show statistical uncertainties. The dashed his- 
togram indicates the signal from a squark with Mg = 150 GeV with arbitrary normalisation. 
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Figure 4: Reconstructed invariant mass distributions in all selection channels with data (points) 
events from 255 pb -1 of e + p collisions compared to SM MC predictions. The method used for 
the reconstruction (M e , M rec , M rec v ) depends on the analysis channel. The error band gives all 
model and experimental systematic uncertainties on the SM prediction (solid histogram) added 
in quadrature. Error bars of data points show statistical uncertainties. The dashed histogram 
indicates the signal from a squark with Mg = 150 GeV with arbitrary normalisation. 
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MSSM scenario with photino-dominated lightest neutralino 




tanp=2, |i=-200 GeV, M 2 =80 GeV, M f = 90 GeV tanp=2, n=-200 GeV, M 2 =80 GeV, M, = 90 GeV 




Figure 5: Exclusion limits at 95% CL on (a) X' nk {k = 1, 2) and on (b) A'^ (j = 1, 2) in a 
phenomenological MSSM with a photino (7) like neutralino (x?)- For comparison, the corre- 
sponding limit from the previous HI analysis [1J is also indicated. Also shown are branching 
ratios to the decay channels considered in this analysis for (c) A' 1U . and (d) values at the 
observed limit. 
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MSSM scenario with zino-dominated lightest neutralino 




tanp=2, [1=200 GeV, M 2 =150 GeV,M, = 90 GeV tanp=2, n=200 GeV, M 2 =150 GeV.M, = 90 GeV 




Figure 6: Exclusion limits at 95% CL on (a) A' llfe (k = 1,2) and on (b) A'^ (j = 1,2) in 
a phenomenological MSSM with a zino (Z) like neutralino (Xi)- Also shown are branching 
ratios to the decay channels considered in this analysis for (c) A' llfc and (d) \' x - x values at the 
observed limit. 
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Figure 7: Exclusion limits (95 % CL) on X' llk for (a) k = 1, 2 and (b) k = 3 as a function of 
the squark mass from a scan of the MSSM parameter space. The dark filled region indicates 
values of the coupling A' llfc excluded in all investigated scenarios whereas the light filled region 
is excluded only in part of the scenarios. Indirect limits from neutrinoless double beta decay 
experiments ((3(30u) [|7l[T0l and tests of charged current universality (CCU) flUEl are also 
shown. For comparison, the corresponding limits from the previous HI analysis [1] are also 
indicated. 
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Figure 8: Exclusion limits (95 % CL) on for (a) j = 1, 2 and (b) j = 3 as a function of the 
squark mass from a scan of the MSSM parameter space as indicated in the figures. The dark 
filled region indicates values of the coupling AV^ excluded in all investigated scenarios whereas 
the light filled region is excluded only in part of the scenarios. Indirect limits from neutrinoless 
double beta decay experiments (J3(30v) [TTlfTOll and atomic parity violation (APV) llBlfTOl are also 
shown. For comparison, the corresponding limits from the previous HI analysis [1] are also 
indicated. 
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mSUGRA V 11k =0.3, tan(3=2, \i<0, A =0 
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mSUGRA V 11k =0.3, tan|3=6, u<0, A =0 
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Figure 9: Exclusion limits (95 % CL) in the mo, mi/2 plane assuming \' llk = 0.3 for (a) tan (3 = 
2 and (b) tan/3 = 6 for k — 1, 2 (hatched region) and = 3 (light filled region). A curve of 
constant squark mass is illustrated for m(d) = 280 GeV. Also indicated are constraints obtained 
by the L3 experiment at LEP flU and the D0 experiment at the Tevatron [6]. The dark filled 
region labelled as "not allowed" indicates where no REWSB solution is possible or where the 
LSP is a sfermion. 
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Figure 10: Exclusion limits (95% CL) in the m ,m 1 / 2 plane assuming A'^ = 0.3 for (a) 
tan/3 = 2 and (b) tan/3 = 6 for j = 1,2 (hatched region) and j — 3 (light filled region). 
Curves of constant squark mass are illustrated for m(u) = 275 GeV and m(t) = 270 GeV. Also 
indicated are constraints obtained by the L3 experiment at LEP [4] and the D0 experiment at 
the Tevatron 0. The dark filled region labelled as "not allowed" indicates where no REWSB 
solution is possible or where the LSP is a sfermion. 
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mSUGRA M = m = m 1/2 , u<0, A =0 




tan (3 

Figure 1 1 : Exclusion limits for M = m = myz in mSUGRA as function of tan (3. Shown are 
the 95 % CL exclusion domains for the model parameters from the production of squarks of first 
and second generation (u, c and d, s) and of third generation (t, b) assuming a value of A' = 0.3 
for the respective coupling. The area below the curves is excluded. The dark filled region region 
labelled as "not allowed" indicates where no REWSB solution is possible or where the LSP is 
a sfermion. 
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